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ABSTRACT 

Using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope, 
we have obtained mid-infrared photometry of the least massive known brown 
dwarf in the Chamaeleon I star-forming region. For this young brown dwarf, 
OTS 44, we have constructed a spectral energy distribution (SED) from 0.8- 
8 fim by combining the measurements at 3.6, 4.5, 5.8, and 8.0 /zm from IRAC 
with ground-based photometry at J, J, H, and K. The resulting SED for OTS 44 
exhibits significant excess emission longward of 3 fim relative to the SED expected 
from the photosphere of the brown dwarf. We have successfully modeled the 
source of this excess emission in terms of an irradiated viscous accretion disk 
with M < 10~^° M0yr~^. With a spectral type of M9.5 and a mass of ~ 
15 Mjup, OTS 44 is now the coolest and least massive brown dwarf observed 
to have a circumstellar disk. These measurements demonstrate that disks exist 
around brown dwarfs even down to the deuterium burning mass limit and the 
approximate upper mass limit of extrasolar planetary companions. 

Subject headings: accretion disks - planetary systems: protoplanetary disks - 
stars: formation — stars: low-mass, brown dwarfs — stars: pre-main sequence 



1. Introduction 

Measuring the lowest mass at which brown dwarfs harbor circumstellar disks is crucial 
for understanding the formation mechanism of brown dwarfs and for determining whether 
planets can form around low-mass brown dwarfs. The natural laboratories for performing 
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this measurement are the nearest regions of star formation (r ~ 0.5-3 Myr, d — 150-300 pc), 
particularly those that have been surveyed extensively for brown dwarfs (e.g., Luhman 1999; 
Comeron et al. 2000; Martin et al. 2001; Briceiio et al. 2002). For some of the young brown 
dwarfs discovered in these surveys^ , indirect evidence of circumstellar disks has been obtained 
through detections of accretion. The most readily observed signature of accretion in a young 
brown dwarf is a broad Ha emission line profile, which has been observed through high- 
resolution spectroscopy in GM Tau, IC 348-382, and IC 348-415 at a spectral type of M6.5 
(White & Basri 2003; MuzeroUe et al. 2003) and in IC348-355 and 2MASS 1207-3932 at M8 
(Jayawardhana et al. 2003b; Mohanty et al. 2003). Very large equivalent widths of Ha in 
low-resolution spectra also have been suggestive of accretion in KPNO-Tau 12 and S Ori 55 
at M9 (Luhman et al. 2003a; Zapatero Osorio et al. 2002) and in S Ori 71 at LO (Barrado 
y Navascues ct al. 2002). High-resolution data has confirmed the presence of accretion in 
KPNO-Tau 12, as well as in several other young sources at M6.25-M8.5 (MuzeroUe et al. 
2005; Mohanty et al. 2005). 

Direct evidence of circumstellar disks around brown dwarfs has been provided by detec- 
tions of infrared (IR) emission above that expected from stellar photospheres alone. Modest 
excess emission in the K and L bands has been observed for several young objects at M6-M8 
and for a few as late as M8.5 (Luhman 1999, 2004; Lada et al. 2000, 2004; Muench et al. 
2001; Liu et al. 2003; Jayawardhana et al. 2003a). At longer, mid-IR wavelengths, the In- 
frared Space Observatory (ISO) detected excesses for CFHT 4 at M7 (Pascucci et al. 2003), 
Cha Ha 1 at M7.75 (Persi et al. 2000; Comeron et al. 2000; Natta & Testi 2001b), and 
GY141 at M8.5 (Comeron et al. 1998). Photometry at similar wavelengths was recently ob- 
tained from the ground with Gemini Observatory for CFHT 4 (Apai et al. 2004), Cha Ha 1, 
and 2MASS 1207-3932 (Stcrzik ct al. 2004). Detections of circumstellar material around 
young brown dwarfs have been extended to millimeter wavelengths for CFHT 4, as well as 
for IC 348-613 at M8.25 (Klein ct al. 2003). Additional detections of accretion and disks 
have been reported for objects in Ophiuchus with IR spectral types of M6 and later (Testi et 
al. 2002; Natta et al. 2002, 2004; Mohanty et al. 2004). However, because those types were 
derived from water absorption bands by using field M dwarfs as the standards, they likely 
have systematic errors such that the temperatures, and thus the masses, are underestimated 
(Luhman et al. 2003b, 2004). 

To continue to search for circumstellar disks around brown dwarfs at later types and 
lower masses, we have obtained mid-IR photometry with the Infrared Array Camera (IRAC; 
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Fazio et al. 2004) on the Spitzer Space Telescope for source 44 from Oasa et al. (1999) (here- 
after OTS 44), which is the coolest and least massive known brown dwarf in the Chamaeleon I 
star-forming region (M9.5, M ~ 15 Mjup; Luhman et al. 2004). In this letter, we describe 
these observations, construct a spectral energy distribution (SED) for OTS 44, measure its 
mid-IR excess emission, and compare this excess to our model predictions for emission from 
a circumstellar disk. 



2. Observations 

As a part of the Guaranteed Time Observations of the IRAC instrument team, we 
obtained images at 3.6, 4.5, 5.8, and 8.0 fim with IRAC on the Spitzer Space Telescope of 
the northern cluster in the Chamaeleon I star- forming region on 2004 July 4 (UT). The plate 
scale and field of view of IRAC are 1'.'2 and 5'2 x 5'2, respectively. The camera produces 
images with FWHM= l'.'6-l'/9 from 3.6 to 8.0 //m (Fazio et al. 2004). The cluster was 
mapped with a 6 x 7 mosaic of pointings separated by 280" and aligned with the array axes. 
At each cell in the map, images were obtained in the 12 s high dynamic range mode, which 
provided one 0.4 s exposure and one 10.4 s exposure. The map was performed twice with 
offsets of several arcseconds between the two iterations. The resulting maps had centers 
of a = 11^^09"^26% 6 = -76°36'26" (J2000) for 3.6 and 5.8 /im and at a = llHO'^'W, 
S = -76°30'20" (J2000) for 4.5 and 8.0 /im, dimensions of 33' x 29', and position angles 
of 28° for the long axes. The images from the Spitzer Science Center pipeline (version 
SI 0.5.0) were combined into one mosaic at each of the four bands using custom IDL software 
developed by Robert Gutermuth. 

Aperture photometry for OTS 44 was extracted with the task PHOT under the IRAF 
package APPHOT using a radius of two pixels (2'.'4) for the aperture and inner and outer 
radii of two and six pixels for the sky annulus. In the next section, we will compare the 
IRAC measurements of OTS 44 to those of KPNO-Tau 4, a young brown dwarf in Taurus 
with the same spectral type as OTS 44. Therefore, we also measured KPNO-Tau 4 from 
images obtained by Hartmann et al. (2005) with the same aperture sizes applied to OTS 44. 
Because the 8.0 /xm images of KPNO-Tau 4 were contaminated by residual images from a 
prior observation of a bright star, the photometry in this band is uncertain and therefore is 
not considered in this work. For an aperture radius of 10 pixels and a sky annulus extending 
from 10 to 20 pixels, we adopted zero point magnitudes (ZP) of 19.601, 18.942, 16.882, 
17.395 in the 3.6, 4.5, 5.8 and 8 /xm bands, where M = -2.5log{DN/sec) + ZP. We then 
applied aperture corrections of 0.210, 0.228, 0.349, and 0.499 mag to the photometry of 
OTS 44 and KPNO-Tau 4. The resulting IRAC photometry for OTS 44 and KPNO-Tau 4 
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are listed in Table 1. The photometric uncertainties for both objects are 0.02 mag at 3.6 and 
4.5 /im and 0.04 mag at 5.8 and 8.0 /xm, with an additional 10% uncertainty in the IRAC 
calibration. Tabic 1 includes optical and near-IR photometry for OTS 44 and KPNO-Tau 4. 
The data at / are from K. Luhman (in preparation) and Briceho et al. (2002) and the near-IR 
measurements are from Oasa et al. (1999) and the 2MASS Point Source Catalog for OTS 44 
and KPNO-Tau 4, respectively. 

3. Analysis 

Before discussing OTS 44, we first examine the SED of the comparison source, KPNO- 
Tau 4. Figure 1 shows the SED of KPNO-Tau 4 derived from the photometry in Tabic 1 
using a distance modulus of 5.76 for Taurus (Wichmann et al. 1998) and the appropriate 
zero-magnitude fluxes (Bessell 1979; Bessell & Brett 1988; Fazio et al. 2004). No correction 
for extinction was applied to these data because of the negligible extinction toward this object 
(Briceiio et al. 2002). As shown in Figure 1, the mid-IR SED of KPNO-Tau 4 is consistent 
with a blackbody at the effective temperature of the central object (Teff = 2300 K; Luhman 
et al. 2004)^, and thus exhibits no evidence of excess emission from circumstellar material. 
Therefore, in the following analysis, we adopt the SED of KPNO-Tau 4 to represent the 
intrinsic photosphere of young objects at the spectral type of OTS 44. 

We now examine the SED of OTS 44, which is shown in Figure 1 after an extinction 
correction of Aj = 0.3 (Luhman ct al. 2004) and adopting a distance modulus of 6.13 for 
Chamaeleon I (Whittet et al. 1997; Wichmann et al. 1998; Bcrtout et al. 1999). In comparison 
to KPNO-Tau 4, OTS 44 exhibits excess flux in the IRAC bands that becomes stronger with 
longer wavelengths. We have compared this excess emission with the predictions of a model of 
an irradiated accretion disk. For these calculations, we adopted a stellar mass = 0.015 Mq 
and radius i?* = 0.23 R© (Luhman et al. 2004), uniform grain size distribution everywhere 
given by the standard power law n{a) ~ a~^"^ with minimum and maximum grain sizes 
flmin — 0.005 /im and a^ax — 0.25 /im (Mathis, Rumpl, & Nordsieck 1977), uniform mass 
accretion rate, and an inner radius = 3 i?* (given by the silicate sublimation radius). 
The disk was heated by local viscous dissipation and stellar irradiation. The emission from 
a wall at the dust destruction radius was included. Detailed descriptions of the method used 
to calculate the disk structure and emergent intensity are given by D'Alessio et al. (1998), 
D'Alessio et al. (1999), and D'Alessio, Calvet, & Hartmann (2001). 



^The mid-IR SED of an object at the age and temperature of KPNO-Tau 4 is predicted to agree with 
that of a blackbody (Allard et al. 2000). 
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As shown in Figure 1, we were able to fit the observed SED with a disk model seen pole- 
on (i — 0) with an accretion rate of M ~ lO~^°M0yr~^ and an inner radius of 3 i?*. We 
were driven to low inclinations because we found that any significant contribution from the 
inner wall, which radiates approximately as a blackbody at 1400 K, resulted in an SED that 
was inconsistent with the IRAC observations. Our estimate of the accretion rate is higher 
than values derived for other young brown dwarfs from modeling of Ha emission profiles 
(MuzeroUe et al. 2003, 2005). Within the context of our model and the adopted inclination, 
we required substantial viscous dissipation due to accretion as well as irradiation from the 
central star to achieve high enough temperatures at the inner edge of the disk. Our model 
assumes that the inner disk wall is perfectly vertical (perpendicular to the disk plane) . If the 
disk wall is actually somewhat tilted or curved, as was adopted to model the older T Tauri 
star TW Hya (Calvct ct al. 2002), also seen pole-on, it is possible that the intermediate 
range of temperatures at the inner disk edge resulting from such structure would enable us 
to fit the observations with lower M. Further exploration of this issue is needed. 

The accretion himinosity for the adopted values of the stellar mass and radius and the 
disk mass accretion rate is 0.0002 L0, which corresponds to 15% of the stellar luminosity, 
L^, = 0.0013 Lq (Luhman et al. 2004). If matter from the disk is loaded onto the star through 
a magnetospheric flow, as the observed line proflles in higher mass brown dwarfs and very low 
mass stars seem to indicate (MuzeroUe et al. 2003, 2005; White & Basri 2003; Jayawardhana 
et al. 2003b; Mohanty et al. 2003, 2005), then an accretion shock should form at the stellar 
surface, where most of the accretion luminosity would be emitted (Calvet & GuUbring 1998; 
GuUbring et al. 2000). This emission would veil the photospheric fines. To estimate the 
expected degree of veiling, we have calculated the emission from accretion shocks on the 
stellar surface following the procedures of Calvet & GuUbring (1998). Figure 1 shows the 
shock emission from two models of accretion columns with values of the energy flux in the 
column = 10^^ erg cm^^ s^^ and T = 10^^ erg cm~^ s~^, which arc representative of those 
estimated for higher mass T Tauri stars (Calvet & GuUbring 1998) and correspond to flUing 
factors on the stellar surface of 0.0016 and 0.00016, respectively. As shown in Figure 1, veifing 
from the shock is predicted to be significant at optical wavelengths but negligible longward 
of 1 /im. Optical spectroscopy can test the model prediction of a high accretion rate; as 
discussed above, modifying the structure of the inner disk edge could result in acceptable 
fits with lower accretion rates and thus lower accretion shock luminosities. 

The IRAG data do not constrain the maximum grain size. Disk models with different 
grain maximum sizes have the same SED for A < 9 /im because the fiuxes arise from the 
optically thick and flat inner disk regions. Thus, observations at longer wavelengths are 
necessary for constraining the dust properties in the OTS 44 disk. Similarly, it is not possible 
to constrain the disk mass with the available data; models with a different viscosity parameter 
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a and disk radius (which imphes different disk masses for a given mass accretion rate) have 
the same SED in the observed spectral range. 

4. Discussion 

Through observations with IRAC aboard the Spitzer Space Telescope, we have shown 
that the least massive known member of Chamaeleon I, OTS 44, exhibits strong mid-IR 
excess emission, which we have reproduced with a model of an accretion disk. With a spec- 
tral type of M9.5 and a mass of ~ 15 Mjup, OTS 44 is now the coolest and least massive 
brown dwarf observed to have a circumstellar disk. The presence of a disk around OTS 44 
demonstrates that the formation of free-floating bodies via disks extends down to the deu- 
terium burning mass limit (M ^ 14 Mjup; Burrows et al. 1997), and raises the possibility of 
planet formation around objects that themselves have planetary masses. Indeed, a candidate 
planetary companion has already been discovered near the earlier and more massive young 
brown dwarf 2MASS 1207-3932 (M8, M ~ 30 Mj^p; Gizis 2002) by Chauvin et al. (2004). 

The relative ease of the detection of the excess emission from OTS 44 (Tint = 20.8 s) 
demonstrates the feasibility of searching for disks around young brown dwarfs at even lower 
masses with IRAC (M < 10 Mjup). In addition, photometry at longer wavelengths with the 
Multiband Imaging Photometer for Spitzer and spectroscopy with the Infrared Spectrograph 
on Spitzer of OTS 44 and brown dwarfs like it should provide detailed constraints on the 
physical properties of disks around low-mass brown dwarfs, such as their mineralogy and 
geometry (e.g., Apai et al. 2004; Furlan et al. 2005). 
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Table 1. Photometry for OTS 44 and KPNO-Tau 4 



ID 


I- J 


J-H 


H-Ks 


Ks 


[3.6] 


[3.6] - [4.5] 


[4.5] - [5.8] 


[5.8] - [8.0] 


OTS 44 


~ 4.7 


1.01 


0.79 


14.61 


13.67 


0.41 


0.53 


0.70 


KPNO-Tau 4 


3.73 


0.97 


0.74 


13.28 


12.51 


0.10 


0.12 
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Fig. 1.— SEDs of the young brown dwarfs KPNO-Tau 4 and OTS 44 (points), (a) The 
slope of the mid-IR SED of KPNO-Tau 4 is reproduced by a blackbody at the effective 
temperature of its stellar photosphere (Teg = 2300 K, solid line), and thus exhibits no excess 
emission from circumstellar material, (b) To estimate the SED of the stellar photosphere 
of OTS 44, the measured SED of KPNO-Tau 4 is combined with the blackbody SED at 
A > 6 iim and scaled to the if-band flux of OTS 44 {short dashed line). The excess flux 
above this photospheric SED is modeled in terms of emission from a circumstellar accretion 
disk (dotted line). The sum of this model disk SED and the photospheric SED agrees fairly 
well with the data (solid line). The values of the model parameters are M — 10~^° yr"-*^, 
an inclination angle to the line of sight i = 0° (i.e., pole-on), a — 0.01, and Rwaii — 3 i?*. 
We also include the SEDs predicted for the accretion shock for values of the energy flux that 
are representative of higher mass T Tauri stars, T — 10^^ and 10^^ergcm~^s~^ (long dashed 
lines) . 



